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Abstract 

A variational and perturbative treatment is provided for a family of generalized 
spiked harmonic oscillator Hamiltonians H = — ^ + Bx 2 + ^ + , where B > 0, A > 
, and a and A denote two real positive parameters. The method makes use of 
the function space spanned by the solutions \n> of Schrodinger's equation for the 
potential V(x) = Bx 2 + . Compact closed-form expressions are obtained for the 
matrix elements <m\H\n>, and a first-order perturbation series is derived for the wave 
function. The results are given in terms of generalized hypergeometric functions. It is 
proved that the series for the wave function is absolutely convergent for a < 2. 
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1. Introduction 

The term "supersingular potential" was introduced by Detwiler and Klauder [1] to 
describe the situation where every matrix element of a perturbation with respect to 
the unperturbed state is infinite. An important class of problems with this feature is 
the family of spiked harmonic oscillators with Hamiltonian H given by Harrell [2], 
where 

d A 

H = -—r+x 2 ^ < x < oo, a,A>0, (1.1) 

dx z x a 



and a > 3. This Hamiltonian has had a variety of applications in atomic, molecular, 
nuclear and particle physics, where it provides a simple model in which the potential 
has a repulsive core \x~ a . Since the early investigations of Harrell, spiked harmonic- 
oscillator potentials have become the subject of an intensive study lasting over two 
decades [3-16]. Most of this work is concerned with approximations for the energy 
of the problem in one spatial dimension, because the Klauder phenomenon [17-19] 
associated with the Hamiltonian (1.1) fails to occur in higher dimensions. Meanwhile, 
no definite results have been obtained so far on a perturbation expansion for the 
wave function associated with the Hamiltonian (1.1), although Aguilera-Navarro and 
Guardiola [4] have reported that they met serious convergence difficulties in the case 
a = 2 with expressions for the first-order correction to the wave function. 

In earlier work [12-14] we have pointed out the advantages of basing our per- 
turbation analysis on an exactly soluble model which itself has a singular potential 
term [12]. We examined a family of generalized spiked harmonic-oscillator Hamilto- 
nians 

H = H a + XV = + Bx 2 + 4 + 4r> B > 0, A>0, (1.2) 

ax 2 - x A x a 

defined on the one- dimensional space (0 < x < oo) with eigenfunctions satisfying 
Dirichlet boundary conditions, that is to say, with wave functions vanishing at the 
boundaries. The singular basis, consisting of the set of exact solutions of H , serves 
as a better starting point for the analysis of the singular problem than the Hermite 
functions used in most earlier work on the Hamiltonian (1.1). In this present paper, 
we derive compact expressions for the matrix elements developed earlier. Moreover, 
we prove convergence of the first order correction to the wave function, using standard 
Rayleigh-Schrodinger perturbation theory for the case a < 2 . Specifically, our main 
results can be summarized as follows: 

1. The matrix elements of the operator x~ a , with respect to the exact solutions of 
the Gol'dman and Krivchenkov Hamiltonian [20-21] H , are given by the following 
explicit expressions 



x -» - (-l)n+rn (|)n ^7 ~ f ) Kl)n{l)m , a a ^ _ n _ —• 1) (13) 

x mn -{ lj a ^ r(7) y n!m , 3 ^2( m,7 2 ,1 2 ,7,1 n 2 ,l), (1.6) 

where 7 = 1 + |vT+4A, valid for all values of the parameters 7 and a such that 
a < 27 . These expressions are also very convenient for direct computations. In 
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particular, these expressions, drawing their importantance from their simplicity, 
are used in studies of so-called 'vestiges' [22] of potentials, left after the cofficient 
A 10. 

2. The first-order perturbation correction of the wave function of Hamiltonian (1.2), 
for a = 2 , is 

4'\x) = ^± T J^^-h e -^ [ln{ ^) _ mi (i.4) 

where ip(l) is the digamma function, also known as the logarithmic derivative of 
the r -function, defined by ip(x) = j^lnT(x) . We note that the coefficient of the 
exponential term is a logarithmic epression in the variable x, instead of being a 
polynomial. This fact explains the odd result that the perturbation expansion fails 
to be regular. In the case of a < 2 , we demonstrate that the first order correction 
of the wave function is 



oo 

x f e VB(c+iy) {c + i y) -7 ( i_^_)3 J r 2 (i,i,i + " ;2 ,2;l--^)dy. 
J c + iy 2 c + iy 



(1.5) 



The functions 3 F 2 are special cases of the generalized hypergeometric function defined 
by 

p 

oo n (atih zk 

pF q (a 1 ,a 2 , . . .,a p ;(3 1 ,p 2 , ■ ■ -,/3 q ;z) = ^ Tf' ^' 6 ^ 

i=i 

where p and q are non-negative integers and (j = l,2,...,q) cannot be a non- 
positive integer. The expression (a) k denotes the Pochhammer symbol 

r(aj + k) 

(a,i) = 1, {a i ) k =a l {a l + l)...(a t + k-l)= ) , k = 1,2, 

1 («i) 

which we also give in terms of the gamma function T(a) . The function defined by 
(1.6) includes as special cases the hypergeometric and confluent hypergeometric func- 
tions 2-F1 and iFi respectively. If the series does not terminate (i.e., none of an, 
i = 1,2, . . . ,p, is a negative integer), then the series converges or diverges according 
as \z\ < 1 or \z\ > 1 . For z = 1 on the other hand, the series is convergent, provided 

q p 

£ ft - £ on > 0. 

3=1 i=l 

2. The matrix elements 

Landau and Lifshitz [21] exhibited a closed form description of the exact solutions 
of the one-dimensional Schrodinger equation (in units % = 2m = 1 ) 

+ (Bx 2 + 4)^n(a;) = E n ip n (x), x e [0,oo), (2.1) 
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where the ip n (x) satisfy the Dirichlet boundary condition ip n (0) = 0. The wave func- 
tions have the form 



Mx) ^ ] n>=(-ir^ ^^ x^h-^\F 1 (-n, 7 ,VBx% (2.2) 
with exact eigenenergies 

E n = 2VB(2n + j), n = 0,l,2,..., (2.3) 

for 7 = 1 + + AA . Herein, 1F1 stands for the confluent hypergeometric function 
defined in (1.6). We have further introduced the alternating coefficient terms (-1)™ in 
the definition of ip n {x) in order to guarantee a smooth transition, through the identity 



2' ' 2x (2n+l)! J 



lFl (-n,-,x 2 ) = V J 77T I ,„ H 2n+1 (x), (2.4) 



to the solutions of the harmonic oscillator problem A = . 

Using the integral representation of the confluent hypergeometric function, Hall 
et al proved that the matrix elements of the generalized spiked harmonic oscillator 

for a < 27 are 



- a _, 1)m+nB ± r( 7 + m) r( 7 -f) ^ fc /m\ r(fc + 7 -f)r(f-fc + n) 
X ™~ { 1} B Vn!m!r(7 + n) I\ 7 ) ^ ' [k J r(fc + 7 )r(f - fc) ' ^ 

These expressions can be readily obtained and further simplified by means of the 
series representation of the confluent hypergeometric function (2.4). Indeed, using 
(2.2) and (2.4), it immediately follows from 



that 



By resorting to the integral representation of gamma function, we get for -f +^+k+l > 
that 

l rr rr v-v- (-m) k (-n)i B%~1 a 



< > = ^ E E ^f^ r( -i + ^ + * + 1) 

-f + o 

(7)« « 



(2.7) 



i W f-*sfi:^ ),r< " f+7+t+i) 



fc=0 



(-TTt)fc 1 

(7)fc fc!' 
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On the other hand, by use of the definition of the Pochhammer symbols and the series 
representation of the hypergeometric function (1.6), the finite sum inside the bracket 
collapses to 

§ — — = § — — r( -2+^) 

= + 7 + fc) 2 F 1 (-n, -| + 7 + fc; 7 ; 1) (2.8) 

after we invoke Vandermonde's theorem [24]. Consequently, we arrive at the matrix 
elements 

<m\x- a \n> = -TT B* ~* V T(-- + 7 + k) (f ~ k)n ( ~ m)fc 1 

Itt n*-* r(7) ff ( m \ r d + n ~ fc ) r (~f + 1 + fc ) 



2 ( 7 )„^ v ' W r( 7 + A)r(f-fc) 

wherein we have used the identity ( ~^ )fc = (-l) fe (™) . After writing T„ and T m for the 
expressions defined by (2.6), this leads us clearly to the matrix elements (2.5). If we 
now factor out m! from the finite sum appearing in equation (2.9), then we are left 
with 



r(fc + 7 - a/2)r(a/2 -k + n) 
r(fc + 7)r(a/2- fc) 

T(z + 7 - a/2)r(a/2 - z + n) 



Aj ( 1} *I(m-fc)l T(k + 1 )T(a/2-k) ~ 



T(z + l)T(m + 1 - z)T(z + j)T(a/2 - z) ' (2 ' 10) 

where /(z) has at most removable singularities at z = a/2 + n + r for r a positive 
integer such that a/2 + n + r e (-1/2, m + 1/2) . This f(z) is therefore holomorphic 
inside and on a simply closed contour C cutting the real axis at the points -1/2 and 
m + 1/2 . As result thereof, the Euler Summation Formula 

m m 

£/(*0 = + (l/2)[/(0) + /(m)]- fQ 1 (t)f'(t)dt, n 1 (t) = 2^ 8 *" (27rl/t) (2.11) 

o o "=i 

holds [23] on the one hand. On the other hand, we may also apply the concept of the 
total residue [23], in particular 

m 

2 f{k) = 2~wiJ n Ctg{nz) ' /(z)rfZ " S ^ esc{7r ct9 ^ z) ' (2A2) 



fe=0 



where J2 Res c stands for the sum of the residues of the given function at all singular 
points inside the closed contour C . If we let C t and C 2 be the arcs of contour C lying 
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above and below the real axis respectively ( -1/2 is the endpoint of C\ and m + 1/2 
that of C 2 ), then 



-dz, 



m+l/2 

• ! 1/2 L L (2.13) 

T(z + 7 - a/2)T(a/2 -z + n) 



f(z) = e 



T(z + l)r(m + 1 - z)T(z + 7)r(a/2 - z) ' 



The two summation formulas (2.11) and (2.13) are already of great theoretical as well 
as practical value, although the function f(z) of (2.10) is still somewhat cumbersome 
since T(z) is determined in terms of improper integrals. Fortunately we are able to 
derive closed-form formulas for these and also for our summation expression (2.6) by 
the use of the higher order hypergeometric function 3 F 2 . Therewith we obtain the 
exact value of the summation (2.6). 

We further simplify the expression (2.9) by taking note of 

(a-k) n - {1 - a)k{a)n 



(1 - a - n) k 
in order to justify the relation 

r -«_I TT B f -| F (7 - f Kf )n f> (-m) fc ( 7 - |M1 - f ) fc 
-2"^ (7)« ^0 (7)*(l-f "»)**! 

1 a -i r(7 — ^)(^)n , a a a 

- 2 T " r ^"" 7 ( 7 )„ 7 - 2 ' 1 - 2 5 7 ' 1 - 2 - 715 ^ 

After substituting for T n and T m , we finally obtain 



(2.14) 



mn 1 j V n!m! (7)^(7) 3 21 ' 7 2' 2' 7 ' 2 ' j ' 1 j 

which is valid for all values of a and 7 such that 7 - f ^ -fc(fc = 0, 1, . . .) . Herein we 
treat the case a = 2 by a limit, as follows: the limit of the hypergeometric function 
3 F 2 in Eq.(2.15) as a^2 implies 

f ( _ 1)m+ n^ g V(7)n(7) m if n >m, 
[ (_l)^»g^ V(y). if m > n. 

Two important observations follow from these results. First, the matrix elements 
of the Hamiltonian (1.2) are given by ( m,n = 0, 1,2, ..,N - 1 ) 



ff mn =< ro|ff |n >^ 2Vs(2n + 7 )<W + A(-l) m +" J B- • / W» T{l 2 )( 2 ) 



J 1 TD I 



3-F 2 (-m,7- -,1- -57,1 - - -n;l) 



(7)«r( 7 ) (2.17) 
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taken over the N -dimensional subspace spanned by basis (2.2). This is highly suitable 
for systematic calculations of the energy eigenvalues through diagonalization of the 
matrix 

/ H 00 H i . . . Hon-i \ 

flio Hii ••• ffiiv-i (218) 

\-ffjv-io Hn-\i ■■■ Hn-in-i/ 
with the aid of a computer. By increasing the matrix dimension N , we can always 
improve these upper energy bounds. 

Secondly the simple poles occurring at 7- f = can be removed for certain values 
of the parameter A , if A = 7 - f ; thus the matrix elements (2.13) turn out to be 



H mn = 2VB(2n + j)S nm + (-l) m+n B% 



(7)n(7) m r(7-f + l)(f)» 
n!m! (7)^(7) (2.19) 

t? ( a a a 

3 F 2 (-m,7- 2' 1 ~ 2' 7 ' 2 ~ ' ^' 



These reduce, as A = 7 - § tends to zero, to 

H mn = 2^B{2n + 7 )J nm + (-1)™+"-^t\/ ilMlk. ( 2 .20) 

I (7) V n!m! 

These expressions therefore, measure the vestigial effects of the remaining interac- 
tions [22] after the perturbation \x~ a is turned off along the path A = 7 - f — > . 
We thus can identify certain paths for a given a, where the matrix elements (2.13) 
always collapse to 2\ r B(2n+~f)5 nm as A — > . Indeed, if we take for example \/A = 7-f , 
then ^ = i(2VA + a - 2) 2 - \ guarantees the existence of the matrix elements (2.13) 
for the a given; and further, as A tends to zero, the expressions (2.13) approach 
2v / B(2n + j)S nm , as desired. There is nothing special about taking the square root of 
A of course, and it is certainly possible to identify other paths by selecting different 
powers of A for the given a . This comment illustrates the fact, that suitable values 
of the parameter A can lead to complete elimination of the vestigial effects of the 
potential as A | . 

3. Perturbation expansion for a < 7 + 1 

3.1. The Eigenvalue Expansion 

A weak coupling perturbation expansion for a class of spiked harmonic oscillators, 
defined by (1.1) with zero angular momentum, was studied by Auguilera-Navarro et 
al [4]. Their results have been expanded and further improved upon by the present 
authors [13] in order to study the Hamiltonian (1.2). For sake of clarity, we shall 
summarize our results here. The standard perturbation theory to the second order, 
of the ground state energy of the generalized spiked harmonic oscillator Hamiltonian 
(1.2), yields the following expansion 

E = E + A < 0\x~ a \0 > +A 2 f_f + ' • ' ' ( 3 ' L1 ) 

n=l ™ 
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where the E n , n = 0,1,2, are the eigenenergies of the unperturbed Hamiltonian 

H given by (2.3) and < 0\x~ a \n > are expressed by means of (2.15) for n = 0, 1,2 

We have demonstrated that the weak-coupling expansion of the ground state energy 
of the Hamiltonian (1.2) leads to 

r( 7 ) i6 7 r 2 ( 7 ) ( 312 ) 

x 4 -F 3 (l, 1, - + 1, - + 1; 7 + 1,2,2; 1) + . . . , 

wherein the generalized hypergeometric function 4^3 is defined by (1.6). It is clear 
from the convergence condition of 4 F 3 in Eq.(3.1.2), that regular perturbation fails in 
the case of a > 7 + 1 . It is also of interest to observe, that the series 4-F 3 (l, 1, f + 1, f + 
1; 7 + 1, 2, 2; 1) diverges as 7 - § —> .This divergence still holds, even if A = 7 - § — > 0, 
because in this case 4-F3 becomes 3 F 2 (1, 1,7 + l;2, 2; 1) , which diverges for all A > . 
More precisely, the hypergeometric series 3 F 2 (1,1,7 + 1;2,2; 1) diverge for all A>-\. 

In the case of a = 2, 4 F 3 (1, 1, § + 1, § + 1;7 + 1,2,2; 1) reduced to 2 Fi(l, 1; 7 + 1; 1) , 
which possesses a closed form sum given by . Thus the weak-coupling expansion 
for this case becomes 



E(a = 2) = 2VB7 + ^A - ^y^ 2 + • • • , (3-1.3) 

which result coincides with the series expansion of E = VB(2 + + A(A + A)) about 
the point A = . 

3.2. The Eigenfunction Expansion 

Several interesting questions have been raised by Aguilera-Navarro et al concern- 
ing the abnormal behavior of the standard weak coupling perturbation expansion in 
connection with attempts to derive the first-order correction to the wave function. 
Indeed, the first correction to the wave function by means of standard perturbation 
techniques, leads to 

00 

n=l ^° * n 

where V n0 and ip n (x) are given by Eq.(2.11) and (2.2) respectively. We have in this 
case 

To obtain a closed form for the sum appearing in Eq. (3.2.2) for a = 2 , we use the rela- 
tion between the confluent hypergeometric function iFi(-n;7+l; •) and the associated 
Laguerre polynomials L^\-) given by [25] 

1 1[ n ' 7 + i ' > r(n + 7 + l) " U " 
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Thus, the sum in Eq. (3.2.2) becomes 



71=1 



^r(n + 7 r n vv ~~ } (3.2.3) 

= r( 7 )[v>( 7 )-to(v^ 2 )], 



where we invoke formula 17 of the Bateman Project [26], page 214, which originally 
was given by Toscano [27]. Consequently, the first correction to the wave function, 
for the case of a = 2 , takes the form 

^\*) = ^ T -^^^-HH^)-m]e-^. (3.2.4) 

In the case of 7 = 3/2 , or A = 0, and B = 1 , Eq. (3.2.4) reads 

4 1) ( a5 ) = 7r-i a; [/n( !C )-^(|)]e-* a / 2 . 

At this point we note that the coefficient of the exponential term is a logarithmic 
expression in variable x, instead of being a polynomial. This fact explains the odd 
result that the perturbation expansion fails to be regular. The difficulties of summing 
the expression that were reported by Aguilera-Navarro et al [4] (§111 & §VI) arose from 
the existence of the alternating coefficient (-1)" in the matrix elements V 0n ■ However, 
we were able to overcome these difficulties by introducing the alternating coefficient 
(-1)™ in the wave function ip n (x) . It is an important result of the present work, that 
this is justified by the smooth transition of the exact solution of the Hamiltonian 
(2.1) to the exact solution of the known harmonic oscillator problem (A = 0) through 
the relation (2.4). 

We now turn our attention to the evaluation of the sum in Eq. (3.2.2) for a < 2, 
given that 7 > § and B is an arbitrary positive number. This requires a suitable inte- 
gral representation of the confluent hypergeometric function iFi(-n, 7 , \J~Bx 2 ) over an 
appropriate contour, in order to interchange summation with integration and thereby 
readily conclude the absolute convergence of the series just mentioned. We find the 
inverse Laplace transform (integral) representation 

c+zoo 

iFi(o,7,>/B«) = S5( 1 "T)M f e ^H-^(l--)- a dt (3.2.5) 

c — ioo 

under the conditions V~B is nonzero real, J?e( 7 ) > 0,c > 0, \arg{\ - f)| < n, with the 
replacement z = x 2 (x real) to be most advantageous for achieving this end. Further, 
we choose c sufficiently large to guarantee that |1 - f | < n and also \arg(l - |)| < -n 
for z = x 2 . Because we integrate over a path beginning and terminating with c - ioo 
and c + ioo respectively for variable t = c + iy , we must note 

1-7 = 1-^^ = 1-^ + ^ (* = * 2 )> 

t cr + if <r + y z c l + y z 
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for which 

t c z + y z 

provided c is choosen large enough - i. e. x 2 < 2c . Moreover, the condition |arg(l-|)| < 
7r must also be satisfied, which translates into \arg(l - ^-)| < n , and this is clearly 
true, provided x 2 < c . Therefore, c must be choosen so large that ox 2 . For such c 
we shall always have 

&■ + y £ 



^„ 2 — V ^i?. 



Now we further continue the process of evaluating the summation in terms of the 
preceeding inverse Laplace- Transform representation (3.2.5) written for a = -n and 
z = x 2 , namely 

c+oo 

1 F 1 (-n, 7 ,v A Bx 2 )=i?5(i-7)^ J e VBt t - 1{l _^ ) n dt (n = 0,1,2,...), (3.2.6) 

c— ioo 

which subsituted into the summation of Eq. (3.2.2) yields 

_ w C+lOO 

/ a \ 1 00 / a \ 1 /■ - 2 

y ^2 ) = B I(l- 7 ) M -l r( } y 12^ _L / e ^*t-7(l _ f.jr.rff 

^-^ n n! ^— ' n n\ t 

n=l n=l „ 

n— 1 ' 



00 ( cx \ 1 r 2 



c + iy 

(3.2.7) 

The evaluation of this last infinite sum, involving integrations over the interval 
( -00,00), is achieved by examining the summation of the integrands, namely 

°° (SL\ 1 r- T 2 1- °° (-) 1 T 2 

yWh^ _L e >^(c+iw)( c+iy )-7(i £_)» = e V5(«H-<»)( c+il ,)-7y- kin_L(i E_)n (3.2.8) 

z — ' n TV. r.4- 7,11 — ' n, n\ c. -4- rii 



71=1 



c + iy n n! c + 



and demonstrating that it has an Li(-oo, oo) -majorant. Hence, the existence of such 
a majorant shall permit us to interchange summation with integration, as result of 
the Lebesgue Dominated Convergence Theorem [28]. To arrive at such a majorant, 
we continue by noting that 

H, = F(f +n) ^ (| + n)f+"-^ e -f-"V2^(l + (l)) 2 
nn\ nr(f)r(n + l) ^ n r(f)(n + l)"+3e-«- 1 V2T(l + o(l)) 

as consequence of the Stirling formula for large arguments of the gamma function. 
By means of cancelations as well as utilizing the fact that 
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we may easily conclude that 

(f)n = (f +n)^r 1 +"e 1 -f 
n n\ ' r(f)(l + n)i+™ 



(i + (i))= (1 r(fP (1 + o(1)) - 



(3.2.10) 



Herein we utilize the well known sequence (1 + ^)™ defining the e p , expressed asymp- 
totically as 



f +nY 
1 + n 

We therefore have that 



= 1 + 



— — 1 \ ™ 



1 + n 



(f)n 



n! r(f)(n+l) 2 "2 



1 (l + o(l)) as n^oo. 



« [1 + 0(1)], 



(3.2.11) 



(3.2.12) 



which makes possible a comparison to the standard p -series ]T ^ < oo provided 

n=l ™ 

p > 1 . We thus have that 



OO /„\ oo ^ 

< A(a) -o— 3- < oo with A(a) < oo 



(3.2.13) 



^(n + l) 2 -f 

for all a such that 2 - § > 1 or a < 2 . We now return to the majorization of 
summmation (3.2.8) in terms of inequality (3.2.13), which entails 

n=l ' " ,V n=l ' y 



< A(a,\/B, c)|c + m/| 
Therein the constant A(a, Vb, c) is defined as 



A(a,VB,c) = A{a)e cV ^ 



E( n+1 ) 4 



n=l 



(3.2.14) 



(3.2.15) 



and also |l - < 1 was made use of. The most important aspect of inequality 
(3.2.14) is the appearance of the L x (-00,00) -function |c+iy|~ 7 ofvariable y majorizing 
the series 

E ^ |e *«) ((+1>) -, (1 _-^, 1 

n=l 

and this aspect justifies the evaluation of summation (3.2.7) by means of the Lebesgue 
Dominated Convergence Theorem [28]. Thus we specifically have 

Y ilk JL lFl (-n r/ , VB.T 2 ) = B^^EM Y ^1 / e^ c +^( C + zy)^x 
z — ' n n! 2-7ri ^ n n\ J 



n=l 



n=l 



(l-__)«idy 
c + 



00 / Q \ 

1 



(1 



-)' 



OO 

B *(1-7)EM£ / e^c+ij/) 

2^ 2 7 



^ n n! c + ly 
2 



(c + - -J— ) 3 F 2 (1, 1, 1 + 2, 2; 1 - -J--)dy, (3.2.16) 

c + ij/ 2 c + ly 
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which is an effective, straight-forward and precise determination of the summation 

J2 ^ IL ^iF 1 (-n,^,^/B~x 2 ) in terms of integrals of higher order hypergeometric fune- 
ral 

tion. However, by utilizing t = c + iy we reconvert the last expression of relation 
(3.2.16) to the inverse Laplace transform format, namely 

oo (SL \ 

{ n n\ 

n=l 

c+ioo (3.2.17) 

c—ioo 

The particality of this result lies in the fact that we can sum our series (3.2.7) involving 
hypergeometric function 1F1 on any finite interval (-c, c) in terms of either the Fourier 
Integral (3.2.16) or (even better) the inverse Laplace- Transform (3.2.17) involving the 
higher order hypergeometric function 3 F 2 as given by the last integral expressions 
(3.2.16) and (3.2.17) respectively. This is exceedingly more effective than bringing 
in the Euler Summation Formula or the summation formula in terms of the total 

residues of 7r[ctg(nz)]f(z) . 

4. Conclusion 

In this paper we have continued our study of the spiked harmonic-oscillator prob- 
lem by expressing the Hamiltonian H as a perturbation of the singular Gol'dman and 
Krivchenkov Hamiltonian H . The zeroth-order eigenf unctions generated by H form 
a suitable singularity-adapted basis for the appropriate Hilbert space. Our principal 
results are threefold: (1) the derivation of the compact closed form (1.3) for the matrix 
elements (2) the addition of a new expansion (1.4) for the wave function; and 

(3) the proof in §(3.2) that the series for the wave function is convergent for all a < 2. 
The compact matrix-element formula will certainly be useful for the generation of 
explicit eigenvalue estimates in particular cases. Aguillera-Navaro and Guardiola [4] 
have discussed the difficulties inherent in deriving a perturbation series for the wave 
function and they have presented a wish list of results to be sought. We are happy to 
report that in the present work we have now found most of the results asked for in 
this list. 
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